NASA/CR-2000-210114 
ICASE  Report  No.  2000-20 


Optical  Measurement  of  the  Speed  of  Sound  in  Air 
Over  the  Temperature  Range  300-650  K 


Roger  C.  Hart 

ICASE,  Hampton,  Virginia 

R.  Jeffrey  Balia  and  G.C.  Herring 

NASA  Langley  Research  Center,  Hampton,  Virginia 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  VA 


Operated  by  Universities  Space  Research  Association 


National  Aeronautics  and 
Space  Administration 


Langley  Research  Center 
Hampton,  Virginia  23681-2199 


Prepared  for  Langley  Research  Center 
under  Contract  NAS  1-97046 


April  2000 


DTIC  QUALITY  INSPECTED  2 


DISTRIBUTION  STATEIV1ENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


20000509  121 


OPTICAL  MEASUREMENT  OF  THE  SPEED  OF  SOUND  IN  AIR  OVER  THE 
TEMPERATURE  RANGE  300-650  K 

ROGER  C.  HART* ,  R.  JEFFREY  BALLA*,  AND  G.  G.  HERRING* 

Abstract.  Using  laser-induced  thermal  acoustics  (LITA),  the  speed  of  sound  in  room  air  (1  atm) 
is  measured  over  the  temperature  range  300-650  K.  Since  the  LITA  apparatus  maintains  a  fixed  sound 
wavelength  as  temperature  is  varied,  this  temperature  range  simultaneously  corresponds  to  a  sound  frequency 
range  10-15  MHz.  The  data  are  compared  to  a  published  model  and  typically  agree  within  0. 1-0.4%  at  each 
of  21  temperatures. 
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1.  Introduction.  The  speed  of  sound  in  air  near  room  temperature  has  been  intensely  studied  over 
a  range  of  various  parameters  such  as  frequency,  humidity,  and  pressure.  This  body  of  work,  stretching 
back  three  centuries,  has  been  reviewed  by  Wong  [1].  Most  of  this  work  was  conducted  in  the  audio  (20- 
20,000  Hz)  or  hundreds-of-kHz  frequency  range.  Less  frequently,  work  was  carried  out  in  the  MHz  range. 
At  room  temperature,  sound  speed  models  for  air  [1]  and  pure  gases  [2]  agree  with  measurement  to  0.03% 
or  better.  However,  previous  work  on  the  speed  of  sound  in  air  is  limited  to  temperatures  of  ~  300  K  and 
below.  In  this  letter  we  present  measurements  of  the  speed  of  sound  (10-15  MHz)  in  air  at  1  atm  over  the 
300-650  K  temperature  range  and  compare  these  data  with  theoretical  values.  This  letter  is  an  outgrowth 
of  our  earlier  work  [3]  in  developing  remote  diagnostics  for  flow  temperature  in  wind  tunnels. 

We  compare  our  experimental  data  to  a  model  synthesized  by  Zuckerwar  [4].  The  speed  of  sound  is 
written  in  the  standard  Laplacian  form  with  corrections  for  various  gas-phase  phenomena.  Thus, 

(1.1)  Vsound  =  [(' VoRT/M)(l  +  <7i)(l  +  c2)(  1  +  C3)]0-5 

where  70  is  the  ratio  of  specific  heats,  R  is  the  universal  gas  constant,  T  is  the  temperature,  and  M  is 
the  molecular  mass.  C\  is  a  correction  for  the  temperature  variation  of  70,  C2  is  a  correction  for  virial 
(non-ideal  gas)  effects,  and  C3  is  a  combined  correction  for  the  vibrational,  rotational,  and  translational 
(classical  thermal  and  viscous)  relaxations.  For  our  experimental  conditions,  C\,  C2,  and  C3  are  each 
<  0.05.  Additional  details  of  these  corrections  are  described  in  Refs.  [3]  and  [4]. 

2.  Experimental  Method.  The  experimental  method  is  laser-induced  thermal  acoustics  (LITA)  [5]. 
Detailed  descriptions  of  the  LITA  apparatus,  procedure,  and  data  analysis  are  described  in  Ref.  [3].  Sound 
waves  are  produced  by  crossing  two  high-peak-power  pulsed  laser  beams  (532  nm)  at  ~  1  deg,  forming 
interference  fringes  at  the  crossing  inside  an  oven.  Two  counter-propagating  ultrasonic  sound  wave  packets 
are  launched  by  electrostriction.  The  propagation  of  these  two  acoustic  waves  is  monitored  by  probing  them 
with  a  long-duration  laser  pulse  (750  nm)  that  is  Bragg- diffracted  into  a  fast  optical  detector.  This  probe 
signal  decays  as  the  wave  packets  propagate  out  of  the  sample  volume  (defined  by  the  probe  beam)  and 
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Table  3.1 

Summary  of  Model  and  Measurement  for  the  Speed  of  Sound  in  Air  from  300  to  650  K  Using  LITA 


T  N  Nr  H  fLTTA  VIATA  Vmod  Ay  oin  (7AV 


K 

% 

MHz 

296.9 

1432 

15 

34 

10.18 

303.8 

445 

6 

33 

10.21 

313.0 

282 

3 

46 

10.55 

321.1 

294 

3 

29 

10.51 

346.0 

292 

3 

30 

10.92 

374.7 

288 

3 

25 

11.29 

381.3 

282 

3 

27 

11.48 

408.4 

277 

3 

43 

12.28 

415.2 

591 

6 

31 

11.95 

431.7 

299 

3 

45 

12.39 

443.1 

297 

3 

26 

12.31 

467.8 

660 

7 

34 

12.75 

482.6 

497 

6 

32 

12.84 

499.3 

903 

10 

32 

13.17 

520.6 

776 

9 

40 

13.57 

528.4 

496 

6 

34 

13.54 

551.4 

257 

3 

36 

13.87 

570.2 

188 

3 

41 

14.25 

603.1 

274 

3 

43 

14.67 

631.3 

257 

3 

41 

14.99 

653.2 

498 

6 

41 

15.29 

m/s 

m/s 

m/s 

m/s 

m/s 

347.5 

346.2 

1.1 

0.018 

0.51 

350.7 

350.2 

0.73 

0.032 

0.75 

354.2 

355.7 

-1.5 

0.043 

0.081 

360.3 

360.0 

0.31 

0.034 

0.035 

374.6 

373.7 

0.92 

0.036 

0.095 

388.7 

388.9 

-0.20 

0.054 

0.066 

392.8 

392.3 

0.46 

0.042 

0.15 

407.6 

406.3 

1.2 

0.060 

0.25 

409.6 

409.5 

0.30 

0.033 

0.14 

415.8 

417.8 

-2.0 

0.053 

0.23 

423.6 

423.0 

0.61 

0.054 

0.14 

434.7 

434.8 

-0.12 

0.041 

0.81 

440.3 

441.6 

-1.1 

0.05 

0.28 

449.2 

449.3 

-0.03 

0.038 

0.84 

458.0 

458.9 

-1.0 

0.043 

1.34 

461.5 

462.2 

-0.58 

0.053 

0.38 

472.7 

472.3 

0.45 

0.071 

0.13 

481.0 

480.5 

0.5 

0.12 

0.18 

492.7 

494.3 

-1.6 

0.10 

0.31 

504.5 

505.7 

-1.3 

0.10 

0.41 

514.7 

514.6 

0.24 

0.089 

1.13 

as  the  waves  are  attenuated  through  relaxation.  In  principle,  the  speed  of  sound  can  be  determined  from 
the  wavelength  (derived  from  beam  crossing  angle)  and  frequency  of  the  temporal  modulation  of  the  signal. 
This  modulation  arises  from  the  beating  of  two  individual  diffracted  light  beams  (one  from  each  sound  wave 
packet)  that  are  Doppler  shifted  relative  to  each  other. 

Since  we  cannot  measure  our  1-deg  crossing  angle  to  0.1%  accuracy  geometrically,  we  effectively  measure 
it  with  another  LITA  apparatus.  We  simultaneously  generate  a  second  LITA  signal  with  the  same  crossing 
angle  in  a  reference  cell  at  known  temperature  of  ~  295  K.  The  sound  speed  in  the  oven  is  then  determined 
from  the  ratio  of  the  LITA  frequencies  observed  in  the  oven,  fLITA,  and  reference  cell,  /re*^ ,  and  the  known 
sound  speed  Vref  in  the  reference  cell.  Oven  temperatures  were  read  with  a  precision  of  0.2  K  with  a 
thermocouple,  which  had  been  calibrated  by  comparison  with  a  precision  NIST-traceable  Pt  thermometer. 
Allowing  for  small  temperature  gradients  in  the  oven,  we  estimate  an  absolute  uncertainty  in  the  measurement 
of  the  oven  temperature  of  0.33%  (±1  K  at  300  K.  growing  to  ±2  K  at  600  K).  The  small  day-to-day  variations 
in  atmospheric  pressure  produce  negligible  changes  in  the  speed  of  sound,  relative  to  the  accuracy  that  we 
achieve  here. 

3.  Experimental  Results.  The  results  of  our  experiment  are  summarized  in  Table  3.1.  In  what 
follows,  quantities  associated  with  a  single  laser  shot  use  subscript  i,  quantities  averaged  over  a  single  run 
(100  shots  in  ~  2  min)  use  subscript  j,  and  all  quantities  appearing  in  the  table  (without  a  subscript  i  or  j) 


2 


are  averaged  over  a  temperature  bin  of  3--15  runs.  A  total  of  107  one-hundred-shot  runs  were  made,  of  which 
more  than  9500  shots  yielded  usable  data.  The  data  shown  have  been  binned  into  21  temperatures,  the 
width  of  the  bins  varying  from  0  K  (all  runs  at  the  same  temperature)  to  5.0  K.  The  reason  for  consolidating 
the  107  runs  into  21  temperature  bins  is  to  reduce  the  length  of  Table  3.1.  For  each  good  shot  of  each  run, 
the  oven  sound  speed  ViLITA  =  {fi'ITA / was  determined,  and  for  each  run,  an  unweighted  average 
y  lit  A  wag  founcL  Using  the  observed  oven  temperature,  humidity,  and  the  unweighted  average  of  measured 
sound  frequency  fjjITA  for  each  run,  we  calculated  yP10^  from  the  model  (Eq.  (1.1)). 
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Fig.  3.1.  Velocity  difference  (AV  from  Table  3.1)  between  the  measurements  and  the  model  as  a  function  of  temperature. 
Uncertainties  due  to  temperature  uncertainties  and  sound  velocity  uncertainties  (am  from  Table  3.1)  were  added  in  quadrature 
to  obtain  the  total  la  uncertainties ,  indicated  by  the  error  bars. 

In  Table  3.1,  Nr  is  the  number  of  runs  and  N  is  the  total  number  of  shots  used  in  each  bin.  The 
average  temperature  T,  relative  humidity  if,  and  model  velocity  Vmod  for  each  bin  were  calculated  as  the 
average  of  per  run  values  weighted  by  the  number  of  good  shots  in  each  run,  e.g.  T  =  'E(NjTj)/1ENj . 
Average  values  of  fLITA  and  yLITA  were  weighted  with  the  per  run  standard  deviations  of  the  mean  <jj. 
The  quantity  erm  was  calculated  from  the  per  run  standard  deviations  of  the  mean  aj  weighted  by  Nj 
according  to  am  =  [ E(Nj(Nj  -  1  -  1))]  /  ;  thus,  om  is  the  standard  deviation  of  the  mean 

of  the  distribution  if  all  the  runs  in  a  bin  are  normalized  to  have  the  same  mean  value.  The  average  of  the 
differences  A Vj  —  y^ITA  —  y™od ?  weighted  by  per  run  standard  deviations,  is  denoted  by  AV,  while  the 
standard  deviation  of  the  AV},  again  weighted  by  per  run  standard  deviations,  is  denoted  by  <jav- 

Comparison  of  the  am  and  ctav  allows  evaluation  of  the  short  (<  10  min)  and  long-term  (>  1  day) 
repeatability  of  the  method.  All  bins  containing  three  runs  comprise  data  taken  over  a  ten  minute  period; 
here  ctav  is  only  a  few  times  am  (average  value  —  2.7  for  Nr  =  3).  Bins  containing  more  than  three 
runs  contain  data  from  more  than  one  day  and  show  substantially  higher  values  of  a  a  v/&m  (~  30  for  the 
Nr  =  15  bin).  Clearly  our  implementation  of  the  reference  cell  technique  does  not  provide  long-term  stability 
comparable  to  the  short-term  precision.  The  uncertainty  due  to  temperature  and  the  statistical  uncertainty 
crm  are  comparable,  and  adding  them  in  quadrature  results  in  the  total  uncertainties  (1  o)  shown  as  error 
bars  in  Fig.  3.1.  Our  measurements  are  in  complete  agreement  with  the  model,  and  significantly  smaller 
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uncertainties  will  be  necessary  to  find  deviations  from  the  model. 

4.  Conclusion.  Our  results  lead  us  to  believe  that  this  noninvasive  method  for  generating  and  detecting 
free-field  ultrasonic  waves  is  promising  for  future  work  in  gas-phase  acoustics.  The  statistics  of  our  sound 
speed  measurements  can  be  improved  by  averaging  more  laser  shots.  In  addition,  improvements  in  future 
LITA  apparatus  and  ovens  are  possible.  Thus,  we  anticipate  that  LITA  can  be  used  for  characterization 
of  the  speed  of  sound  for  a  variety  of  gases,  frequencies,  and  temperatures  with  improved  uncertainties 
compared  to  those  we  have  achieved  here. 

We  gratefully  thank  A.  Zuckerwar  for  helpful  discussions. 
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